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Abstract 

A lattice theory of scalar bosons in the fundamental representation of the gauge 
group SU (Nc) and of the global symmetry group SU (Nf) is shown to induce a standard 
gauge theory only at large Nf. The system is in a deconfined phase at strong scalar 
self-coupling and any finite Nf . The requirement of convergence of the effective gauge 
action imposes a lower limit on the scalar mass. 
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1 Introduction 



There is ample numerical evidence, but still no analytic proof for confinement in nonabelian 
gauge theories. Lattice simulations of pure gauge theories clearly show that large Wilson 
loops follow an area law, (W) ~ e~^^, which is usually taken as the signal for confinement. 
Besides the scale parameter (or coupling constant) a nonabelian gauge theory contains only 
one parameter, the number of colors, N^. There were numerous attempts to explore simplifi- 
cations which might occur at the large Nc limit and allow for analytic solutions. Confinement 
in lattice gauge theories can be proven at strong coupling but in the relevant, weak coupling, 
scaling, region no analytic approach has succeeded. 

Kazakov and Migdal [|I| conjectured that a gauge theory induced by a self adjoint scalar 
field might be exactly solvable in the Nc ^ oo limit. An induced gauge theory has a simpler 
structure because its bare lagrangian does not contain the kinetic and self-interaction terms 
of the gauge fields, familiar from QCD. The hope was that after integrating out the scalar 
bosons one ends up with an effective gauge theory, which is in the same universality class 
as a standard lattice gauge theory. The pure gauge terms of the lagrangian are supposed to 
be generated by scalars loops. The theory would become solvable due to the existence of a 
"master" scalar field. It has been pointed out, however that the lagrangian of this model 
possesses an extra gauge symmetry, which makes the expectation value of the standard 
Wilson loop vanish identically 0. Such a property shows that unlike in QCD the quark- 
antiquark potential is infinite. Quark-antiquark systems are pointlike. 

Various solutions have been offered to the problem of extra gauge symmetry. "Filled 
Wilson loops" (combinations of loops in which every link variable and its conjugate appear 
the same number of times) as the measure of confinement and the perturbative (in variable 
l/Nc) introduction of pure gauge terms have been proposed 0. 

A way to eliminate the extra gauge symmetry is to introduce scalar fields in the funda- 
mental representation The lattice lagrangian for a scalar field with compact gauge 
interaction is linear in the gauge field, U, and does not have the extra gauge symmetry 
of the lagrangian for adjoint representation fields. Simulations have been performed with 
interacting fields in mixed representations @] and the large limit has been investigated P|. 

Our aim is to investigate an induced gauge theory with fundamental scalars in the Nf 
dimensional fundamental representation of a global flavor group SU{Nf). Rather then taking 
the limit Nc oo, we will try to use analytic approximations at finite Nc- This approach 
is similar to that of Hasenfratz and Hasenfratz , who investigated a gauge theory induced 
by fermions in the fundamental representation, where the lagrangian is also linear in the 
compact gauge field. One of their conclusion, based on perturbative calculations, was that 
both in the continuum and lattice theories a minimal number of fermions, Nf > llNc/2, is 
needed to achieve asymptotic freedom in the induced gauge theory. They also pointed out 
that in the limit of large fermion mass one needs to take the limit of Nf — oo to be able to 
obtain a conventional lattice gauge theory. 

Using fundamental scalars instead of fermions makes the problem more complex, due to 
the self-interaction of scalars, which is required by renormalization theory. As we will see 
later the scalar self coupling has a significant effect on the phase of the induced gauge theory. 
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A perturbative one loop level analysis shows that the induced theory is asymptotically 
free only if the number of flavors, Nf, satisfies the inequality Nj > 22Nc. This condition is 
very similar to that of Hasenfratz and Hasenfratz in fermion induced gauge theories . 
The main purpose of the present paper is to investigate confinement, or rather the conditions 
under which the induced gauge theory is in the same universality class with a standard gauge 
theory with a Wilson plaquette action. 



2 Lattice gauge theory with fundamental scalars 

In the absence of the pure gauge (plaquette) term the lattice lagrangian density for 
scalars in the fundamental representation of the gauge group, SU{Nc) and of the global 
flavor symmetry group, SU {Nf has the following form: 



2D 



L = m^(Pl{x)(P\x) + Ai + A2 ' + ^ + e„)f/n(x)0'(x). (1) 

n=l 



where the summation is over all 2D unit lattice vectors (positive for n < D and negative 
for D < 77, < 2D) originating at x. Note that f/„(x) = Ul_^_4^{x + in), where n = 1,...,D. 
Summation over gauge indices has been omitted, i and j represent flavor indices. We will 
use a rescaled version of this Lagrangian, redefining the scalar fields (f){x) as m(f){x) — > (j){x). 
We will also introduce rescaled couplings Xijvn? Aj. This rescaling will amplify the role 
of the large m limit that will be repeatedly used below. Note that one can only expect a 
confined induced theory in the momentum range below the scalar mass, m serves as a cutoff 
in the theory . We will also see that the condition of convergence of the effective gauge 
lagrangian imposes a lower bound onm/n, as well. 

We will use a transformation to make the lagrangian quadratic in the scalar field. In- 
troducing sources ji{x) and j2{x) where the first is a gauge scalar, while the second is a 
hermitian gauge matrix, we can write the generating functional as 

Z = exp\ -Ai r . / X2 - ^2Tr ^ . \ / d4){x)dUi{x) 

exp I - ^ [0^(x)(ji + j2 + l)0(x) + —Y^cf)\x + e,)[/,(x)0(x)] |. (2) 

Note that if Nc > Nf then the tensor = (f)f(f)^^, where a and (3 are gauge indices, is 
degenerate, consequently j2/3 also satisfies some constraints. In what follows, we will assume 
that Nf > Nc. 

The integral over the scalar field, being Gaussian, can be performed to give 

Z = expj-AiE-^-A2Tr^-^| ( dU{x)d-'' , (3) 
i X ^3Ax)^ x^32[x)^]J 

where we dropped the subscript / from Nf. 
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Now the determinant, d, can be rewritten using Schwinger's proper time formalism as 



d = Det 



K 



j{x)+jo{x) + 1 + —U{x) 



-N 



exp < A^Tr 



-T\j(x)+jo{x)+l+^Uix)] 

T 



(4) 



Here f/ is a step operator, connecting nearest neighbor lattice points. U satisfies the following 
constraint: {x\U\x') = {x'\W\x). Note that the singularity at T = is an irrelevant constant, 
independent of U and jj. After expanding the exponential in (^) and transforming the 
operator products into functional integrals one obtains the following representation for 
the generating functional 



Tr exp < — / dr, 



Jl{x^'\r.))+J2{x'''{r^)) + ^U{x,{r,)) 



K 



dx^'\n) 

(5) 



where Tj is the length of the ith path. Each of the path integrals is path ordered in the 
SU{Nc) space. The n paths are given in parametric representation, x^^\t,i), satisfying the 
condition xW(Ti) = a;«(0). 

Each term of the sum over n in (|^) contains n path integrals over closed loops. We will 
evaluate these functional integrals below. Notice that while the operators ji{x) are diagonal 
in coordinate space, operator U generates discrete steps on the lattice. For a finite size loop 
there can only be a finite number of discrete steps. Only a finite power of the combination 
[ti^ /m?)U will contribute. Expanding one of the functional integrals of (|^) we obtain 



dx(r)Tr exp 



dr 



Ji{x{t)) + J2{x{t)) + —U{x{t)) 



EETr 



exp {-ro [ji(xo) + 32{xo)]} \{ Ui-u exp {-n [ji{xi) + j2{xi)]} 



i=l 



(6) 



where summation over £k denotes all possible lattice paths of length k lattice units. Sum- 
mation over Ti, the length of the T-interval between two subsequent operators, Ui-u and 
Uii+i, can be turned into integrals over Ti, /HiLidrj. The variables r^ satisfy the relation 
J2i=o'Ti = T. Introducing a multiplier 1 = / dro5(ro + ti + ... + t^ — T) we can see that 
the integrand depends only on ro + r^. Integrating over ro introduces a multiplier ro + r^. 
After the substitution tq + Tk t^ and adding all cyclic permutations of the operators Ui 
we obtain for (H) 



]J f/i_iiexp {-ri [ji{xi) + j2{xi)]} 



i=l 



(7) 



The multiplier -1 of jvn^ was omitted because k, the number of links of a closed loop is 
always even. 
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Notice now that exp{— r[ji +72]} are translation operators for the functional derivatives 
with respect to the sources ji and 72 • After translating the differentiations and setting ji = 
in (^) we obtain extra terms in the effective action. These extra terms represent contact path- 
path interactions. The U matrices are path ordered and traced along n different paths. The 
contact interaction proportional to Ai is scalar, but the other one, proportional to A2 is 
proportional to operator Pjj that crosses the SU{Nc) indices of path i and path j at the 
point of intersection. The integration over T can be performed using the delta-functions in 
(0) to result in the following expression for the generating function 



n=0 



n 

k=l 



(k) 



n 



k,l=l 



E 
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(8) 



The jth vertex of the fcth loop is denoted by xf\ rj'^^ denotes the "time" (fraction of T) 
spent at that vertex. The partition function is a sum over the products of n closed loops of 

In other words, 
is a loop-loop 



perimeter r^. Cr^. symbol is used to identify specific loops of perimeter r^. 



Cr^ defines the shape and location of the loop. The last multiplier in 
interaction term generated by the scalar self-interaction. 
If the scalar couplings vanish then the last multiplier of 



is 1 and the integrations over 



Tj and the summation over n can be performed to result in the exponentiated form 



/rff/.(x)exp|iV;^(^^^ ETrnf^,,+i|. 



(9) 



This is exactly the form of the partition function one would obtain in a much simpler way, 
namely by expanding the logarithm of the boson determinant in a power series of [t^^ /m?)U 
and summing over cyclic permutations of the product Tr 11^=1 Uj. 

A necessary condition for the convergence of the expansion can easily be obtained. Take 
an arbitrary closed loop. Now split the loop open at any of its vertices, Xq, and add an 
"appendix" consisting of a set of k [/-matrices on links, {xi, Xj+i), i = 0, /c — 1, ending up 
at vertex Xk- Close the loop by adding all the k links to end up at xq. In other 

words, the new loop, obtained after this surgery, is doubled up on the segment between 
Xq and Xk- Due to pairwise cancellations the corresponding contribution of U matrices is 
{k^ /rn?Y^Uxo,xi---Uxk-ik^Xk,Xk-i---Uxi,xQ = (K^/m^)^*^. Consequently one obtains a contribu- 
tion of the same form as that of the unsplit loop. Clearly, the expansion of the action in (^) 
converges only if the sum over all the possible additions of the above kind converges. Since 
the doubled up path is a random walk of length /c, the number of these walks is {2D)'' the 



condition for convergence is m > \'2Dk . This condition gets modified if Aj 7^ 0. Then no 
exact condition can be obtained, because self-intersections of the random walks between xq 
and Xk alter the weight. Neglecting the effect of self- intersections of the doubled up line the 
condition for convergence get modified as 



m 



(10) 
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where 

/■°° 1 1 /1\ 

Ux) = / dridrae-"^-"^-"^"^^ = -e'Ei - . (11) 
JO X \xj 

The origin of the notation I2 and further properties of this function will be discussed later. 
Note that every projection operator P generates a multiplier Nc by creating one more closed 
gauge index loop along the doubled up path. The function hix) satisfies the condition 
-^2(2;) < 1- -^2(0) = 1. The coefficients of the weak coupling expansion are all finite but the 
expansion, having a zero radius of convergence, is only asymptotic. I2{x) — (l/x)logx if 
X is large. The physical reason for the decreasing nature of F[x) is clear: the scalar self- 
interaction and with that the path-path interaction is also repulsive, making the excursion 
from xq to Xk energetically less and less favorable. In the limit of A, — » 00 doubled up loops 
are forbidden. In general, the larger A is, the better the convergence of the expansion. 

A few more comments about (H) are in order. The contribution k = I and i = j renormal- 
izes m and can be omitted. This can be seen in the representation of the partition function 
without rescaling the scalar field, m(j) — »■ 0. The contributions i ^ j or/and k ^ I oi the 
contact term prevent us from resumming the series over n. The summation can only be 
performed using a linked cluster expansion. One expects the induced gauge theory to be 
QCD-like only at momentum scales below m, when the only role of the scalars to generate 
a gauge coupling, because free bosons cannot be produced 

Let us now examine the effect of path-path interactions more closely. Our aim is to 
generate an effective action using a linked cluster decomposition. We will call a certain set of 
n lattice loops a diagram. Two loops are connected if they have at least one common vertex. 
A subdiagram is connected if there is a continuous path formed from links of its loops between 
every pair of its vertices. The contribution of connected diagrams to the effective action can 
be calculated in the following manner. Suppose that there are Sm vertices in the diagram 
where m loops intersect. A point where m loops intersect (including self intersections) 
contributes to the diagram by 

TTL f 771 I 

/™ = / n [dTiC-^^] exp \ - ^(Ao + Pij\i)nT, \ (12) 

Then the contribution of a connected diagram consisting of n loops of perimeter r^, 
k = 1, ...,n to the effective action equals 

S = N-{^] n Tr n CiW^^^'-''^^ (13) 



k=l j=l 



where 



Wii,-,k) ^ Yl^^y^ (14) 



m 



is the weight factor appearing due to self and mutual intersections of k connected loops. 
The linked cluster expansion has the following form 



k=l ri,r2,.-.,r-fe V"' / Conn[£ri ■■■•Crfc] «=1 



X(i.2,...,fc)^ (15) 
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where Conn[Crj^Cr2---jCrf.] means that summation is over connected diagrams only. The 
coefficients X*^^'^' "'^^ satisfy the following recursion relation 



_ ^1 J2 6{k-mi-2m2- ...-kmk)f[ 



k 1 rj^(l,2„...,j) 

(16) 



rrti,...,mfe i=l 

We can solve these recursion relations iteratively to arrive at the following expression for 
= 1, 2, and 3: 

X« = W^^\ (17) 

X(i'2) = iy(i'2)_iy(i)iy(2)^ (18) 

X(1.2,3) _ ^^(1,2,3) ^ 2Vr(^)l^(2)|^(3) _ ^(1)^(2,3) _ ^(2)^(1,3) _ |^(3)y^(l,3)^ ^^g) 

Of course, the last three terms of X^^'^'^^ lead to identical terms in 5*. 

At infinite coupling one only has self avoiding walks. Then the overlap functions W^^''"'^'^ = 
for A; > 1 and W^^^ = 1. The recursion relation ([IB]) can be solved exactly to give the 
following weights 

XiK:k) ^ [ ^) (20) 

rC 

The worthiness of the cluster expansion depends on the value of the scalar couplings. In 
fact, for large couplings the cluster coefficients, X^^'\ tend to a finite constant. Contributions 
for diagrams having multiple overlaps decrease very rapidly. The overlap function - 
has the following asymptotic behavior: 

^(1,...,.) ^ aogf\ (21) 

A* 

where t = J2m=2^Sm/'2- Thus, it is advantageous to use expansion with respect to 
If the coupling is small, however, then the cluster expansion is more appropriate, because 
— I _ 0(A), while the asymptotic behavior of cluster coefficients is 

^(l,...,r) _ (22) 



3 Strong scalar self coupling 

At infinitely strong scalar self coupling the magnitude of the scalar field is fixed. It was shown 
0] that a pure scalar self interacting theory on a lattice at infinite coupling is equivalent to a 
self avoiding random walk. In our system, as shown by (H), the Wilson loops in the partition 
function avoid each other and themselves (note that trivial self interaction of lattice points 
has been previously removed). If no intersection is allowed the integration over r/'^'* can be 
trivially performed. One obtains for the partition function 

z = /d^.(.)E^EE-E £^ E'E'■••E'^Tr^f^f^ (23) 

n=0 ri r2 r„ \ / £ £ fc=l i=l 
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where the primes over the summation signs imply summation over non-intersecting loops 
only. If intersections are not allowed then each gauge variable U^'^^ appears only once under 
the integral and the integral vanishes for all n except for n = 0. Then the partition function 
itself is 1. The expectation value of a Wilson loop, Hi Ui is 

R / 9 \ S 



{I[U^)=N 



(24) 



This is clearly a perimeter law, implying that there is no confinement in the A — oo limit. 
The system is in a deconfined phase. 

The system remains in a deconfined phase at a large enough value of the scalar coupling 
constants. If the coupling is strong we can estimate the expectation value of the loop by 
employing an expansion in 1/A, using the asymptotic behavior of the intersection weights, 
Vl^(A). They behave roughly as l/A*"/^ for a r-fold intersection, as shown by (|2T|) . Then the 
leading contribution to the expectation value of a Wilson loop of area A and perimeter S is 
still (p4D. Contributions of more then 1 loop, having r intersections are cut off by the rth 
power of /2K^/m^ ~ (logA/A)(K^/m^). In particular, the contribution coming from single 
plaquette loops only is 



(Tr n U^) 



K 



(25) 



clearly negligible compared to the main contribution, (|21|) , if N is kept finite at large values 
of A. Also, corrections to this large A expansion converge rapidly. The conclusion is that 
the system is in a deconfined state at sufficiently strong scalar coupling. If is allowed to 
increase as well then the situation may change, as discussed in the following section. 



4 Large N limit 

Let us consider first A, = 0, for simplicity. If Aj 0, than the coefficients = 1 and the 
integrand of the partition function exponentiates. Loops of all sizes contribute as shown by 
(y). It is easy to estimate the contribution of loops of different sizes in the effective action 
to the expectation value of the Wilson loop if H?/m? is small. A loop coinciding with the 
Wilson loop contributes by the same as (P^. A plaquettes, where A is the area of the loop 
contribute by 

S / ^2 \ 4A 



(nC/.l-iV^^-j , (26) 

a negligible contribution, unless A^ is also large. If A^ is small then at vanishing self coupling 
of the bosons the Wilson loop follows a perimeter law. The system is in a deconfined state. 

As shown by (^ we recover the standard Wilson action if we take the simultaneous limit 
of m? — > oo, A^ — ^ oo, such that 1/ g"^ = N {kappa? /m?Y tends to a finite value If tv? 

is not infinite, just large, then one has the area law only for sufficiently small loops: 

A log(mV/t^) 

5 log^^ • ^''^ 
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where A is the area of the loop. We will return to the case of finite N and K^/m^ in the next 
section. 

In the rest of this section we will investigate a theory obtained from (H) by taking the 
limits ^ oo, ^ oo at fixed g'^ = m^/NK,'^. We will allow for arbitrary Aj. Then 
only contributions containing plaquettes survive. The contribution of a connected diagram, 
consisting of two plaquettes, to the effective action is 

S, = [(I,y^ - I] (^1^ TtUp,TtUp„ (28) 

where S2 is the number of points where the two plaquettes intersect (1,2, or 4). The multiplier 
{hY^ — 1 vanishes at Aj = and goes to -1 if A^ oo. 

Let us examine the effect of small couplings, Aj. The cluster functions X*^^''"'^) decrease 
as A**"^. Then obviously two plaquette clusters dominate. One can substitute a pair of 
plaquettes with a cluster of two plaquettes in the covering of Wilson loop approximately 
A different ways. There are four different relevant 2 plaquette clusters, two of which have 
one common point, S2 = 1, and two have two common points, S2 = 2. Altogether one has 
obtains a correction term A[2l2 + 2J| — 4]. Exponentiating the correction term leads to a 
similar correction term to the area law at vanishing A. At small A the correction term is just 
—6XA. Since the correction term is negative, the confining linear potential gets stronger at 
increasing values of A. Clusters of larger size will not change this conclusion, because the 
cluster functions are proportional to higher powers of A and their sum also exponentiates to 
an area law. 

In principle, the cluster expansion could diverge at some finite value of A and fixed g. 
Then the system could go into a deconfined phase. It is easy to see, however, that this is 
not the case. Examine the theory at large A. Then again it is not convenient to use the 
cumulant expansion. Though the gauge coupling might be small, in effect one can use a 
strong coupling expansion, because the overlap of plaquettes is costly, each one introducing 
an multiplier of I2 ~ log A/ A. Then the A plaquette contribution dominates, just as at large 
gauge coupling, giving the following estimate for the expectation value of a Wilson loop: 

{Trnc<>~(^)'. (29) 

For a minimal covering an I4 overlap functions appears at every internal lattice point of the 
Wilson loop. This is so because there are four adjacent plaquettes at every point. Since 
I4 ~ 1/A^, ( p9l) provides a good approximation if g'^X^ » 1. The conclusion is that the 
system is in the confined phase at any fixed g provided A >> 1/g. Since we saw earlier that 
it is in a confined phase at small A as well, it follows that the limit ^ 00, g'^ fixed leads to 
a confining theory at all values of A. When A — 00 the string constant tends to infinity and 
we approach a superconfined theory, similar to the one induced by adjoint representation 
bosons. 
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5 Conclusions 



We have investigated the phases of an SU{Nc) lattice gauge theory induced by complex 
bosons. At the perturbative level such a theory is asymptotically free if the the number of 
flavors, N, satisfies > 22Nc 0. The convergence of the induced theory was shown to 
require a minimum value for the scalar mass. Also, on physical grounds it is expected that a 
sensible induced gauge theory is obtain at energy scales below the scalar mass, which plays 
the role of a cutoff. For these reasons the scalar mass was assumed to be large. 

The self-coupling. A, of scalars, inducing the theory plays an important role. For large 
A and finite but large scalar mass (i.e., fixed "magnitude" for the scalars) the system was 
shown to be in a deconfined state. An obvious limit where a confining theory might be 
obtained is N,m? — > cxd, N{K,'^/m'^)^ fixed. We show that both at small and large A one 
obtains a confining theory (provided the standard lattice gauge theory is confining). Then 
it follows that the limits — > oo and A — oo do not commute. 

It is of considerable interest to determine the phase of the system at finite but large m, 
finite A and N. We have a rather roundabout argument that, at least at some values of the 
parameters, the induced gauge theory has a confined phase. The effective gauge coupling 
(the coefficient of the plaquette action) is ~ Nn'^/m'^, but there are terms in the effective 
action corresponding loops larger then plaquettes. Therefore, it is not at all obvious whether 
the theory confines. However, if m is sufficiently large, then at fixed g coefficients of terms 
corresponding to larger loops decrease rapidly. Imagine now that instead of an induced gauge 
theory we consider a real gauge system with gauge coupling g. The only difference between 
the effective actions of the standard and induced theories is a change in the coefficient of the 
plaquette term, N{k'^ /rn^Y ~^ N{k'^ /rn^Y + l/fi'^- If parameters are such that this change is 
sufficiently small, then no change is expected in the phase of the system because the ratio of 
the plaquette term and terms corresponding to larger loops hardly changes. A real scalar- 
gauge boson system on a lattice is expected to have both confined and deconfined phases and 
at sufficiently large bare scalar mass it will be in the confined phase. Ultimately, the phase 
of the system at finite values of the parameters can only be determined by simulations. 
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